Decoherence of molecular wave packets in an anharmonic potential 
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The time evolution of anharmonic molecular wave packets is investigated under the influence of 
the environment consisting of harmonic oscillators. These oscillators represent photon or phonon 
modes and assumed to be in thermal equilibrium. Our model explicitly incorporates the fact that in 
the case of a nonequidistant spectrum the rates of the environment induced transitions are different 
for each transition. The nonunitary time evolution is visualized by the aid of the Wigner function 
related to the vibrational state of the molecule. The time scale of decoherence is much shorter than 
that of dissipation, and gives rise to states which are mixtures of localized states along the phase 
space orbit of the corresponding classical particle. This behavior is to a large extent independent of 
the coupling strength, the temperature of the environment and also of the initial state. 

PACS numbers: 3.65.Yz, 33.80.-b 



I. INTRODUCTION 

The correspondence between classical and quantum 
dynamics of anharmonic systems has gained significant 
attention in the past few years |], ^, [| . A short laser 
pulse impinging on an atom or a molecule excites a su- 
perposition of several stationary states, and the resulting 
wave packet follows the orbit of the corresponding clas- 
sical particle in the initial stage of the time evolution. 

However, the nonequidistant nature of the involved en- 
ergy spectra causes peculiar quantum effects, broadening 
of the initially well localized wave packets, revivals and 
partial revivals |l], [|, ||, 0, [|, @]. Partial revivals are in 
close connection with the formation of Schrodinger-cat 
states, which, in this context, are coherent superpositions 
of two spatially separated, well localized wave packets 
. Phase space description of vibrational Schrodinger- 
cat state formation using animated Wigner functions can 
be found in |||. However, these highly nonclassical states 
are expected to be particularly sensitive to decoherence 
H . The aim of this paper is to analyze the process of de- 
coherence for the spontaneously formed Schrodinger-cat 
states in the anharmonic potential. 

We consider the decoherence model which relies on the 
interaction of the quantum system (S) with its environ- 
ment (.E) that has many degrees of freedom. The dynam- 
ics of the environment together with the quantum system 
under investigation is unitary, i.e., the density operator 
Pse of the coupled S+E system always represents a pure 
state. Starting from an initially uncorrelated density op- 
erator, the interaction builds up S-E entanglement and 
the reduced density operator 



p s = Tr B (pse) 



(1) 



where Tr^ means a trace over the environmental degrees 
of freedom, turns into a mixture. Assuming negligibly 



short relaxation times in the environment, one obtains 
a Markovian master equation which is a useful tool for 
the dynamical investigation of the environment induced 
decoherence. It allows for the determination of pointer 
stat es llCfl , that is, states that are favored by decoherence 
pj |12[, and also for the calculation of the characteristic 
time of decoherence for different initial states |l2|, |l3| . 
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In the following we introduce a master equation that 
takes into account the fact that in a general anharmonic 
system the relaxation rate of each energy eigenstate is 
different. This master equation is applied to the case of 
wave packet motion in the Morse potential that is often 
used to describe a vibrating diatomic molecule. Con- 
sidering the phase-space description of decoherence, we 
show how the phase portrait of the system reflects the 
damping of revivals in the expectation values of the po- 
sition and momentum operators due to the effect of the 
environment. We also calculate and plot the time evolu- 
tion of the Wigner function corresponding to the reduced 
density operator of the Morse system. The Wigner func- 
tion picture visualizes the fact that although our mas- 
ter equation reduces to the amplitude damping equation 
Jli] , [lj , fL6[ [l7t in the harmonic limit, the anharmonic 
effects lead to a decoherence scheme which is similar to 
the phase relaxation Jljj of the harmonic oscillator (HO). 
It is found that the reduced density operator that arises 
due to the decoherence can be identified with a mixture 
of states that are well-localized in the phase space and 
equally distributed along the orbit of the corresponding 
classical particle. We illustrate the generality of this de- 
coherence scheme by presenting the time evolution of an 
energy eigenstate as well. We also calculate the decoher- 
ence time for various initial wave packets. We show that 
decoherence is faster for wave packets that correspond 
to a classical particle with a phase space orbit of larger 
diameter. 
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II. DESCRIPTION OF THE MODEL 

The total Hamiltonian of an anharmonic system and 
its environment has the form 

H S +E = H S + H E + V. (2) 
The self-Hamiltonian of the system is written as: 

H S = Y, E n\<t>n){K\, (3) 
n 

where the spectrum {E n } is assumed to be nondegenerate 
and discrete, but not necessarily equidistant. Eq denotes 
the ground state energy of the system, and E m > E n 
whenever m > n. The environment is represented by a 
set of harmonic oscillators 

H E =£)fej fc (4a fc + l/2). (4) 

k 

We assume the following interaction Hamiltonian 

V = hX* 9ka k + hX^g k a< k , (5) 

fe k 

where X is an appropriate system operator that trans- 
forms each eigenstate of H$ into a superposition of dif- 
ferent eigenstates corresponding to lower energy values. 
This is the application of the rotating wave approxima- 
tion (RWA) to an anharmonic, multilevel system. The 
operator X* is the Hermitian conjugate of X, and for the 
sake of simplicity the coupling constants g k were taken 
to be real. In the case of a vibrating diatomic molecule 
He describes the modes of the free electromagnetic field, 
while the terms X and X* in the interaction Hamiltonian 
are related to the molecular dipole moment operator. 

We assume that the environment is in thermal equilib- 
rium at a given temperature T, and the environmental 
oscillators have continuous distribution with frequency 
dependent density of states, D(uj). Starting from the 
von Neumann equation for the total density operator, 
standard techniques jl8], |l9[ [H] lead to the following 
Markovian master equation in the Schrodinger picture 

j t Ps = ~ [H S , Ps] ~ (X^X ePs + XXlp s 

- Xl Ps X-X ePs X^ +h.c). (6) 

The matrix elements of the operators appearing in the 
nonunitary terms are given by 

(<Pm\X a \<f> n ) = n(cj) m \X\(j) n ) D(uj nm )g 2 (uj nm )n{uj nm ), 

x {n{uo nm ) + 1) , (7) 

where oj nm = \E m - E n \/%, n(u) = (exp^ - de- 
notes the average number of quanta in the corresponding 
mode of the environment, and the subscripts e and a refer 
to emission and absorption, respectively. 



As we can see, the matrix elements (0) of the operators 
that induce the transitions depend on the Bohr frequency 
of the involved transition, which is a genuine anharmonic 
feature. In the special case of the HO, when H$ has 
equidistant spectrum, and X is identified with the usual 
annihilation operator a, both X a and X e are proportional 
to X = a, and Eq. dq) reduces to the amplitude damping 
master equation Jl4 ]l5|, at a finite temperature. 

In certain cases one can further simplify Eq. (Jfy . When 
the environment induced relaxation rates are much lower 
than the relevant Bohr frequencies, the system Hamilto- 
nian induces oscillations that are very fast even on the 
time scale of decoherence and vanish on the average. Ig- 
noring these fast oscillations we arrive at the interaction 
picture master equation 

-Tj.{<Pi\ps\<t>j) = $i,j^2lik{(t>k\Ps\(t>k) - r C ji{<f>i\Ps\<l>j), 

k^i 

(8) 

that has already been obtained in Refs. (22[ in order 
to treat the spontaneous emission of a multilevel atom. 
In Eq. (H), -fik denotes a relaxation rate, that is the prob- 
ability of the \<pk) \4>i) transition per unit time, while 
^ = l l 2 Y,khik+1jk), where 

f 2 (cpilXelfajifclXlfa) if Kk, 
l lk = { if i = k, (9) 

[ 2 (&|# o |fo)(&|#|0 Jk ) ifi>*. 

However, due to the elimination of the fast oscillations 
related to H$, Eq. (||) is not suitable for investigating 
the wave packet motion and decoherence simultaneously, 
therefore we propose to use Eq. (|^). On the other hand 
we note that Eq. (||) radically reduces the computational 
costs of calculating the time evolution for long times, 
which might be necessary when the system-environment 
coupling is very weak. 

Supposing that our knowledge is limited to the popu- 
lations P n — {4>n\Ps\4>n), both Eq. (H) and Eq. (||) leads 
to the Pauli type equation 

j/n = fr"* Pfe _ IknPn) ■ (10) 
k 

Requiring the condition of detailed balance p4| in 
Eq. ( |Io| ) leads to the steady-state thermal distribution 
at the temperature of the environment. 

When a diatomic molecule is considered in the envi- 
ronment of the free electromagnetic field, the operators 
X and X^ in Eq. (Q) gain a clear interpretation: in the 
eigenbasis of H$ they are the upper and lower triangular 
parts of the molecular dipole moment operator, We 
will assume that ft is linear p5| , that is, proportional 
to the displacement X of the center of mass of the di- 
atomic system from the equilibrium position. In this case 
g 2 (io)D{uj) oc uj 3 , that is, 

{(t>m\X a \<t>n) = A (cj) m \X\(f) n ) UJ nm 3 n{uj nm ), 

(<pm\X e \4> n ) = A (cj) m \X\(j) n ) u) nm 3 (n{u> nm ) + 1) (11) 
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where matrix elements of X can be calculated using 
the algebraic method summarized in [ p6| , and A = 
irg 2 (ui)D(ui)/Lu 3 is an overall, frequency independent cou- 
pling constant. 

However, in order to get insight into the interplay be- 
tween wave packet motion and decoherence, it is worth 
considering a stronger molecule-environment interaction 
than the electromagnetic field modes can provide. Keep- 
ing the structure of Eqs. (pi]), this can be done by in- 
creasing the value of A. Here we present calculations 
with two different coupling constants, Ai and A 2 which 
are chosen so that at zero temperature W01/701 ~ 10 5 
and 4 x 10 3 for A = Ai and A2, respectively. This model 
allows for the numerical integration of the master equa- 
tion (tw (that provides more details of the dynamics than 
Eq. (p|)) in a time interval that is long enough to identify 
the effects of decoherence. These effects can be summa- 
rized in a decoherence scheme (see Sec. |V|) that has a 
clear physical interpretation, and which is valid also in 
the weak molecule-environment interaction, when (||) is 
more efficient to calculate the time evolution. 

In the following we shall apply our master equation 
(^) to the case of the Morse Hamiltonian [^7| that is of- 
ten used to describe a vibrating diatomic molecule. This 
Hamiltonian has the dimensionless form 

H S = P 2 + (s + l/2) 2 [exp(-2A) -2exp(-X)], (12) 

where the shape parameter s is related to the dissoci- 
ation energy Ddiss, the reduced mass of the molecule 
to, and the range parameter of the potential a via 
1/2. The dimensionless displacement and 



_ \/2mD di , 
ha 



momentum operators obey the canonical commutation 
relation [X, P] = i. 

The Hamiltonian (|2|) sustains [s] + 1 normalizable 
eigenstates (bound states), corresponding to the eigen- 



values E m (s) 



-(s — to) z 



0, 1, ... [s], where [s] 



denotes the largest integer that is smaller than s. The 
continuum above E m —\ a ] corresponds to the dissociated 
molecule with positive energies. For the sake of definite- 
ness we have chosen the NO molecule as our model, where 
s = 54.54. 

The initial wave packets of our analysis will be Morse 
coherent states |xo,po), which are localized on the 
phase space around the point (xq,Po)- The Wigner func- 
tion of a representative Morse coherent state is shown in 
Fig. |3] a). These states can be prepared by an appro- 
priate electromagnetic pulse that drives the vibrational 
state of the molecule starting from the ground state into 
an approximate coherent state. An example can be found 
in [E8|, where the effect of an external sinusoidal field is 
considered. 

Let us note that the construction given in [p9| would 
allow us to use arbitrary initial states, but for our current 
purpose it suffices to consider states \xq,po) with negli- 
gible dissociation probability, i.e., coherent states that 
practically can be expanded in terms of the bound states 
|</>„), n = 0,1,... [s\. 



III. TIME EVOLUTION OF THE 
EXPECTATION VALUES 

Starting from \ip(t = 0)) = \xq,pq = 0) as initial 
states, the qualitative behavior of the expectation value 
(X)(t) — (ip(t)\X\ip(t)) draws the limit of small oscil- 
lations. In the absence of environmental coupling (i.e., 
A = 0), for x < 0.05, (X)(t) (as well as (P)(t)) exhibits 
sinusoidal oscillations. For larger initial displacements 




t/t 

FIG. 1: Phase portrait corresponding to the time evolution 
of the initial state \ip(t = 0)) = |a;o,0), with xo = 0.5. The 
parameters are A = Ai, T = 5 TiuJoi/k, and to is the period 
of the small oscillations in the potential. The initial point 
(X) — 0.5, (P) — together with the starting direction is 
also indicated. 



from the equilibrium position, the anharmonic effects 
become apparent: the amplitude of the oscillations de- 
creases almost to zero, then faster oscillations with small 
amplitude appear but later we re-obtain almost exactly 
(A)(0) (and (f)(0) as well), and the whole process starts 
again [0. Without the influence of the environment the 
main attributes of (X)(t) and (P)(t) can be explained by 
referring to the various Bohr frequencies that determine 
their time dependence: dephasing of these frequencies 
leads to the collapse of the expectation value, and we 
observe revival when they rephase again. 

For the initial state of \tp(t = 0)) = \xq,Q), with xq = 
0.5, the original phase of the eigenstates is restored 0, 

around the full revival time t rev = 110 to, where to 
is the period of the small oscillations in the potential. 
At t/to = 55 and t/to = 27.5 half and quarter revivals 
|l|, H can be observed. Fig. [I] shows the damping of 
the revivals both in (X)(t) and (P)(t) when interaction 
with the environment is turned on. Note that the phase 
portrait of the corresponding classical particle would be 
a helix with monotonically decreasing diameter, revivals 
are of quantum nature. However, Fig. [l] does not provide 
a complete description of the time evolution in the phase 
space, this can be given by using Wigner functions, see 
Sec. 
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IV. DECOHERENCE TIMES 

Our master equation (^) describes decoherence as well 
as dissipation. However, the time scale of these processes 
is generally very different, and one can distinguish the 
stages of the time evolution that are dominated either by 
decoherence or dissipation 12 . In Fig. || an example is 



Trp| 




T, the decoherence time is longer for smaller initial dis- 
placements. Additionally, for fixed A and T the func- 
tion td{xo) can be well approximated by an exponen- 
tial curve td(xo) — td(0) exp(— kxq). E. g., for A = Ai, 
T = 10 TbWoi/k and < xo < 2 the parameters take the 
values td(0) = 93 to and k = 0.97. 

It is known that quarter revivals in an anharmonic 
potential lead to the formation of Schrodinger-cat states, 
i.e., states that are superpositions of two distinct states 
localized in space || as well as in momentum || |30| . On 
the other hand, smaller initial displacements correspond 
to classical phase space orbits with smaller diameter. 
Consequently the quantum interference related to non- 
classical states that are formed during the course of time 
cover a smaller area in the phase-space in this case. This 
means that our result is a manifestation of the general 
feature of decoherence that increasing the "parameter of 
nonclassicality" , which is the diameter of the correspond- 
ing classical orbit in our case, causes faster decoherence 
0. A similar result was found in Jl^] for the case of 
decoherence in a system of two-level atoms |D| 31 . 



FIG. 2: The entropy and the purity of the reduced density 
matrix of the Morse system as a function of time, calculated 
using Eq. (^|). The coupling parameter (see Sec. O) is A = Ai 
and T = 10 hujoi/k. The initial state was \4>(t = Oj) = \xq, 0), 
with :ro = 2.0. 

depicted showing how the method of time scale separa- 
tion works. We have calculated the entropy 

S=-Tv[p s Hps)}, (13) 

as well as the quantity Tr[p|], which measures the purity 
of the reduced density operator. Note that the Tr opera- 
tion without subscript refers to the trace in the system's 
Hilbert space. Decoherence time td is defined as the time 
instant that divides the time axis into two parts where 
the character of the physical process is clearly different. 
Initially both S[t) and Tr[p|(t)] change rapidly but hav- 
ing passed td (emphasized by a vertical line in Fig. |^), 
the moduli of their derivative significantly decrease. Af- 
ter td the entropy and the purity change on the time scale 
which is characteristic of the dissipation of the system's 
energy during the whole process |l2|, pH| . In other words, 
decoherence dominated time evolution turns into dissipa- 
tion dominated dynamics around td- In the next section 
we shall determine the density operators into which the 
process of decoherence drives the system. In connection 
with these results we have verified that the states around 
the decoherence time do not change appreciably in a time 
interval that covers the possible errors in determining td- 
An interesting question is the dependence of the de- 
coherence time on the initial state of the time evolu- 
tion. We calculated td as a function of the initial dis- 
placement for the case of displaced ground states (that 
is, coherent states with zero momentum, |xo,0}) as ini- 
tial states. It was found that for all values of A and 



V. WIGNER FUNCTION DESCRIPTION OF 
THE DECOHERENCE 

In order to visualize the time evolution of the reduced 
density matrix of the Morse system we have chosen the 
Wigner function picture, which represents ps as a func- 
tion over the classical phase space 

W(x,p,t) = ±- f (x-u/2\p s (t)\x+u/2)e iup du. (14) 

This description allows us to investigate the correspon- 
dence between classical and quantum dynamics. 

First we consider the ideal case without environment. 
Then, in the initial stage of the time evolution, the posi- 
tive hill corresponding to the wave packet \xo,Po) follows 
the orbit of the classical particle that has started from 
(xq,Po) at < = 0. However, due to the uncertainty re- 
lation, the Wigner function as a quasiprobability distri- 
bution has a finite width, and this fact combined with 
the form of the Morse potential implies the stretching 
of the Wigner function along the classical orbit in the 
course of time. (See Ref. for similar results with the 
Husimi Q function.) After a certain time the increas- 
ingly broadened wave packet becomes able to interfere 
with itself, and around the quarter revival time one can 
observe two positive hills chasing each other at the oppo- 
site sides of the classical orbit. The strong oscillations of 
W between the hills represent the quantum correlation of 
the constituents of this molecular Schrodinger-cat state 
0. Later on the initial Wigner function is restored al- 
most exactly and Schrodinger-cat state formation starts 
again. Detailed Wigner function description of these pro- 
cesses that are related to the free time evolution can be 
found in ||. 
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a) b) 




FIG. 3: Time evolution of the Wigner function corresponding to the initial state \ip(t = 0)) = |aco, 0), with xq = 0.5. The 
coupling parameter is (see Sec. [o]) A = A2 and T = 0.3 TkJoi/k. The plots a) and b) correspond to time instants t/to = and 27.5, 
while both c) and the contour plot d) are snapshots taken at t/to = 137.5. 



In the case when environmental effects are present, we 
found that decoherence follows a general scheme. A rep- 
resentative series of Wigner functions is shown in Fig. |^. 
The snapshots correspond to the initial state and time 
instants when the first and third Schrodinger-cat state 
formation would occur in the absence of the environ- 
ment. Consequently, the Wigner function in Fig. |^ b) 
corresponds almost to a Schrodinger-cat state, but this 
state is already a mixture. However, there are still neg- 
ative parts of the function in between the positive hills 
centered at X\ = 0.51, p\ = and x-2 = —0.34, p 2 = 0. 
The "ridge" that connects these hills along the classical 
orbit is absent in a pure Schrodinger-cat state. Later on 
this ridge becomes more and more pronounced and at the 
decoherence time we arrive at the positive Wigner func- 
tion of Fig. U c) and d). According to the contour plot 
Fig. H d), the highest values of this function trace out the 
phase space orbit of the corresponding classical particle. 



That is, Pg ec , the reduced density matrix that arises as 
a result of decoherence, can be interpreted as a mixture 
of localized states that are equally distributed along the 
orbit of the corresponding classical phase space orbit. 

It is worth comparing this result with the case of the 
HO, when the master equation (@) reduces to the am- 
plitude damping equation |l4|, 15^ 16, |T^], see Sec. ||. 
It is known that harmonic oscillator coherent states are 
robust against the decoherence described by the ampli- 
tude damping master equation (as well as against the 
Caldeira-Leggett master equation jll)), the initial 
superposition of coherent states turns into the statistical 
mixture of essentially the same states. This is a con- 
sequence of the facts that these states are eigenstates 
of the destruction operator a, and the operators in the 
nonunitary terms of Eq. (^) are proportional to a and a? 
in the harmonic case. None of these statements can be 
transferred to the anharmonic system, where the Morse 





coherent states do not remain localized during the course 
of time, even without environment. Therefore the scheme 
of the decoherence is qualitatively different for the har- 
monic and anharmonic oscillators: Our results in the an- 
harmonic system are similar to the phase relaxation in 
the harmonic case [|l4|, where the energy of the system 
remains unchanged, but the phase information is com- 
pletely destroyed. We note that a similar result was ob- 
tained in Rcf. p3| , where the rotational degrees of free- 
dom were considered as a reservoir for the harmonic vi- 
bration of hot alkaline dimers. 

Our decoherence scheme is universal to a large ex- 
tent. In the investigated domain of the coupling con- 
stants Ai < A < A2 and temperatures ranging from T = 
to T = 15 huiQi/k, it is found to be valid for all initial 
states, not only for coherent states. Fig. [| shows an ex- 
ample when the initial state is not a wave packet, it is the 



fifth bound state, corresponding to E 5 , which is very close 
to (0.5, 0|i?s|0.5,0), so direct comparison with Fig. || is 
possible. As we can see, although the two Wigner func- 
tions are initially obviously very different, they follow 
different routes (that takes different times) to the same 
state: Fig. || c) and Fig. [| c) are practically identical. 
The final plot in Fig. [| indicates how the Wigner func- 
tion represents the long way to thermal equilibrium with 
the environment: the distribution becomes wider and the 
hole in the middle disappears. 

It is expected that the loss of phase information has ob- 
servable consequences. According to the Franck-Condon 
principle, the absorption spectrum of a molecule around 
the frequency corresponding to an electronic transition 
between two electronic surfaces depends on the vibra- 
tional state. The time dependence of the spectrum 
should exhibit the differences between the pure state of 
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an oscillating wave packet and the state pg and the 
thermal state. More sophisticated experimental methods 
based on the detection of fluorescence Q or fluorescence 
intensity fluctuations J35|, surely have the capacity of 
observing the dephasing phenomenon considered in this 
paper. 

VI. CONCLUSIONS 

We investigated the decoherence of wave packets in the 
Morse potential. The decoherence time for various initial 



states was calculated and it was found that the larger 
is the diameter of the phase space orbit described by a 
wave packet, the faster is the decoherence. We obtained 
a general decoherence scheme, which has a clear physical 
interpretation: The reduced density operator that is the 
result of the decoherence is a mixture of states localized 
along the corresponding classical phase space orbit. 

This work was supported by the Hungarian Scientific 
Research Fund (OTKA) under contracts Nos. T32920, 
D38267, and by the Hungarian Ministry of Education 
under contract No. FKFP 099/2001. 
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